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ABSTRACT 

We present the discovery of photometric variations in the carbon-dominated atmosphere (hot DQ) 
white dwarf star SDSS J133710. 19-002643.6. We find evidence for two low-amplitude, harmonically- 
related periodicities at 169 s and 339 s, making it the fastest and smallest amplitude of the four known 
hot DQ variables and the only variable whose spectrum suggests the presence of hydrogen. The star's 
fundamental and harmonic amplitudes are roughly equal, and its pulse shape is similar to the other 
two members of the class with detected harmonics. Like the other variables, it appears relatively 
stable in frequency and amplitude. 
Subject headings: stars: individual (SDSS J133710. 19-002643.6, SDSS J142625. 71+575218.3, SDSS 

J220029.08-074121.5, SDSS J234843.30-094245.3) — stars: oscillations — white 

dwarfs 



1. INTRODUCTION 

The McCook-Sion white dwarf catalog contains more 
than 150 white dwarf stars with molecular or atomic 
carbon lines in their spectra (DQ white dwarfs; 
iMcCook fc SionI [1999 and online updates). These 
lines often dominate DQ spectra, yet they are typi- 
cally modeled by atm ospheres containing mostly helium 
(IDufour et al.1 [20051) . Recentlv. however, calculations 
by IDufour et al.l (j2007[) revealed that a small number 
of the hottest DQ white dwarf stars have spectra that 
are fit best by carbon-dominated model atmospheres. 
These hot DQ white dwarf stars have shown themselves 
to be intriguing in other regards as well. It is not, 
for instance, clear where they came from, though re- 
cent theoretical work suggests a scenario in which they 
descended from the hydrogen-deficient PG 1159 stars 
(jAlthaus et al.|[200l [Corsico et al.|[2009[ ). 

Adding to their intrigue and also promising to shed 
light on their origin and nature, recent observations 
have uncovered variables among the hot DQ white 
dwarf stars. Compelled by the idea that carbon- 
atmosp here white dwarf st a rs mig ht be pulsationally un- 
stable, [Montgomery et all (|2008h observed six hot DQ 
white dwarfs to look for photometric variability and 
discovered significa nt periodic luminos ity variations in 
SDSS J1426-h5752. IBarlow et al.1 ([200l then found vari- 
ability in two more hot DQ wh ite dwarfs (SDSS J2 200- 
0741 and SPSS J234 8-0942). IGreen et al.l ^M^ and 
IDufour et al.l (|2009aD present follow-up observations of 
these three stars. 

We here introduce a fourth hot DQ variable, SDSS 
J133710. 19-002643.6 (hereafter SDSS J1 337-0026), a 
star first identified as a DQ white dwarf by ILiebert et al.l 
([2003i) and one of the nine hot DQ stars discussed by 



Based on observations at the SOAR Telescope, a collabora- 
tion between CPNq-Brazil, NOAO, UNC, and MSU. 

t Visiting astronomer, Cerro Tololo Inter-American Observa- 
tory, National Optical Astronomy Observatory, which is oper- 
ated by the Association of Universities for Research in Astron- 
omy, under contract with the National Science Foundation. 



IDufour et all (|2008afn . SDSS J1337-0026 is the fastest, 
smallest-amplitude variable of the four. We detect two 
harmonically-related frequencies: a 339 s fundamental 
and a 169 s first harmonic of comparable amplitude (^ 
0.3%). Their phase relationship is like those of SDSS 
J1426-H5752 and SDSS J2200-0741 and leads to both a 
deep primary minimum and a secondary minimum at the 
location of the fundamental maximum. 

2. OBSERVATIONS 

We observed SDSS J1337-0026 (g'=18.7, u'-g'= -0.46) 
using two different instruments on the 4.1-m SOAR tele- 
scope on Cerro Pachon in Chile. For the first two sets 
of observations , we used the Goodman Spectrograph 
(jClemens et al.l :2004'). an imaging spectrograph mounted 
at one of the SOAR Nasmyth ports. The spectrograph's 
collimator and camera optics produce an image on a 4k x 
4k Fairchild 486 back-illuminated CCD. The plate scale 
at the detector is 0.15 arcsec pixel" ^. 

During the 2009 June and July observations, the spec- 
trograph was not available, so we used the SOAR Optica l 
Imager fSOI: IWalker et al.l [20031: ISchwarz et all [200l . 
The imager is mounted at a bent-Cassegrain port of the 
telescope and has a mosaic of two e2v 2k x 4k CCDs 
read out through four amplifiers. The optics produce a 
plate scale at the detectors of .0767 arcsec pixel" ^. 

We first obtained usable data on SDSS J1337-0026 on 
2008 July 27 during an engineering night for the spectro- 
graph. We observed a field containing the target with the 
Goodman Spectrograph in imaging mode with the CCD 
readout binned 2x2, yielding 0.3 arcsec pixels. We 
collected 1.6 hrs of usable unfiltered photometry during 
which the average seeing was 2.8 arcsec and became in- 
creasingly unstable. To decrease the read-out time, we 
obtained only a 300 x 650 pixel region of interest (ROI). 
On 2009 April 20, we observed SDSS J1337-0026 for 3.6 
hrs with the same setup as before but through a broad- 
band blue S8612 filter and with a 600 x 280 ROI. The 



^ Listed there with erroneous positive declination: 
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SDSS 



Dunlap, Barlow, Clemens 



Table 1 

Observation Log 



Date 


Start Time 


J- exp 


T 

-■- cycle 


Length 


Airmass 


Instrument 


Filter 


Comparison Stars 


(UTC) 


(UTC) 


(s) 


(s) 


(s) 










2008 Jul 27 


00:00:59.4 


45 


49.5 


5939 


1.34-2.01 


Goodman 


none 


Cl'',C2b 


2009 Apr 20 


05:24:24.9 


40 


44.0 


12878 


1.19-3.11 


Goodman 


S8612 


C1,C2 


2009 Jun 27 


23:57:00.3 


38 


39.9 


13390 


1.16-2.09 


SOI 


SDSS g' 


C1,C2,C3<= 


2009 Jul 23 


23:54:00.4 


25 


27.9 


9947 


1.28-2.92 


SOI 


SDSS g' 


C1,C2,C3 



"SDSS J13:37:13. 17-00:28:33.3 g'=16.9 u'-g'=1.5 
''SDSS J13:37:12. 50-00:26:11. 7 g'=17.0 u'-g'=1.2 
'^SDSS J13:37:18. 44-00:25:58.3 g'=17.1 u'-g'=1.9 

average seeing was 1.3 arcsec, and the second half of the 
data, as the Moon was rising, shows significant periodic 
variations in the sky brightness with a period of around 
720 s, as might result from passing clouds with periodic 
structure. 

The next two observing nights presented us with more 
stable atmospheric conditions. On 2009 June 27, we ob- 
tained 3.7 hrs of data through a Sloan g' filter on SOI. 
The whole chip was read out and binned 6 x 6 to yield 
0.46 arcsec pixels, which oversampled the poor average 
seeing of 2.5 arcsec. We gathered 2.8 hrs of data with 
SOI on 2009 July 23 with a ~ 5 min gap due to a guid- 
ing problem. The 4x4 binning resulted in a plate scale 
of 0.31 arcsec pixel"^, and the average seeing was 1.1 
arcsec. Table [T] summarizes these observations. 

3. REDUCTION AND ANALYSIS 

We used IRAF to subtract the bias level and apply 
flat field corrections to the Goodman data. The unfil- 
tered dome flats used to correct the first data set clearly 
do not correct obvious fringing structure present in the 
observation frames. This fringing introduces errors into 
the photometry for that night, but its consequences are 
abated somewhat by the stability of the fringing pattern 
during the observations and the small, incremental drift 
of only about one pixel during the observing run. 

The data from SOI are written in mosaic format since 
each image consists of data from four different amplifiers 
with distinct bias levels and gains. We used IDL to trim 
the overscan regions from each of the four data sections 
and write each exposure into one single-format FITS file 
(ignoring the gap between the two CCDs). Once this was 
done for all the object, bias, and flat frames, we reduced 
the data using IRAF. The resultant images display resid- 
ual differences in mean bias level among the regions read 
out by the four amplifiers, but since we only care about 
differential photometry, it seemed superfluous to employ 
more manipulative means to remove this artifact. 

We analyzed the data by first using a differential pho- 
tometry program we wrote in IDL that uses t he func- 
tion APER, which is based on DAOPHOT (jStetsonI 
Il987f l. We computed signal-to-noise (S/N) estimates 
with different-sized apertures and used the a perture 
that maximizes the S/N in the light curve (jHowelll 
[1989) to perform aperture photometry using the program 
CCD_HSP, an IRAF routine written by Antonio Kanaan 
that also uses DAOPHOT. 

Th e program WQED (v. 2.0, [Thompson fc Mullallvl 
l2009f ) was used to correct the times from UTC to the 
barycentric Julian ephemeris date and to carry out the 
following steps. To mollify the effects of atmospheric 
transparency variations and extinction, we divided the 



hght curve of SDSS J1337-0026 by the average of the 
light curves of comparison stars (Table [!}, which we 
checked against each other for signs of variability. We 
fit a parabola to this divided light curve to approximate 
residual atmospheric extinction effects. Dividing by this 
fit and subtracting one yields the light curves in terms of 
fractional variation about the mean presented in Figure 

dii. 

Though no signal is apparent to the eye in the pho- 
tometry, the amplitude spectra (Fig. [T|d) produced from 
discrete Fourier transforms of each of the light curves re- 
veal noticeable signals on all four nights near 5900 /iHz 
(169 s). A peak near the harmonically-related frequency 
2950 /iHz (339 s) does not always stand out, but the 
combined 2009 amplitude spectrum shows both to be 
obviously above the noise (Fig. (Tt). 

We assess the significance of these peaks from the com- 
bined 200 9 data using the Lornb-Sca rgle normalized peri- 
odogram (; Lomblll976l : IScargld[l98l comput ed by IDL's 
LNP TEST (a routine based on fasper from lPress et al.l 
I1992D . The largest peak in the power spectrum pro- 
duced from the three nights of 2009 data combined is 
near 5900 /iHz and has a power of ~ 26. If we expect 
from the 2008 data a peak near ^ 5900 /xHz and con- 
sider just this frequency, then the probability calculation 
is straightforward. When the periodogram is normal- 
ized by the sample variance, the distribution of powers 
is described by the regularize d incomplete beta function 
(iSchwa rzenber g-Czernvl [l998f ) from which we find that 
the probability that a peak as large as the one near ^ 
5900 /iHz would occur there by chance is ~ 3xl0~^^. 

If, on the other hand, we want to know the false 
alarm probability, i.e., the odds of a peak so large oc- 
curring by chance at some frequency in the range consid- 
ered (0-11370 /iHz, the Nyquist frequency on the April 
night), then we need to know the number of indepen- 
dent frequencies that serve to increase the probabilistic 
resources and thus increase the odds of finding a large 
peak due to noise. We follow the method laid out in sec- 
tion 3Al_a£dAppendi^xBof[CummingIer^ ()1999[ ) (see 
also lHorne fc Baliuna£l ll986') and perform 10^ bootstrap 
Monte Carlo trials. For each trial we compute the power 
spectrum of a light curve constructed with the same ob- 
servation times as the original but with the flux values 
for each time drawn randomly (with replacement) from 
the original flux values. A fit to the high probability end 
of the resulting distribution of maximum powers indi- 
cates that the number of independent frequencies is ~ 9 
times the number of data points (9x961). This results in 
a false alarm probability of~ 3x10^*, which is a factor 
of 10 smaller than an extrapolation of the Monte Carlo 
results. 
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Figure 1. (a) Diflferential light curves for each night and (b) amplitude spectra resulting from their Fourier transforms. Dotted lines 
indicate the positions of the two detected frequencies, (c) Top; Amplitude spectrum (produced with Period04) from the three nights of 
2009 data combined. Significant peaks are present near the harmonically-related frequencies of 2950 /xHz (0.26%) and 5900 /xHz (0.33%). 
Bottom: Amplitude spectrum of the residuals that result from subtracting the best-fit model (with fi «i 2950.856 /^Hz and f2 = 2fi) from 
the combined light curve. No peaks above 700 /xHz are >4 times the mean noise level of 0.055%. 

Similarly, the probability that a peak as large as the 
one near 2950 /iHz {^ 16) will occur there by chance is 
~ 7xl0~^. The probability of a chance occurrence of a 
peak at least this high somewhere in the frequency range 
is ~ 0.07% according to the analytic calculation, or ^ 
0.15% according to the Monte Carlo results. 

To characterize these variations, we model them as 
a sum of sine waves and determine the best-fit ampli- 
tude, frequency, and phase by non-linear least-squares 
fits to the data, which we performed using both Period0 4 
(jLenz fc Bregeii [2001 and MPFIT d Markwardtl [20091 ) . 
The largest peak in the 2008 July amplitude spectrum 
is at approximately half the 2950 /iHz frequency but is 
not significantly present on the subsequent nights, so we 
include this extra, low-frequency component in the fit for 
that night only and note that fitting without it does not 
yield a significant difference. We list the best-fit param- 
eters and their formal errors in Table [2j 

Both SDSS J1426-h5752 and SDSS J2200-0741 have 
non-sinusoidal light curves because of the presence of a 
fundamental and first harmonic. In both of those stars, 
the harmonic minima coincide with the maxima and min- 
ima of the fundamental. This results in a pulse shape 
with a deepened minimum at the location of the funda- 
mental minimum; relatively large harmonic amplitudes 
(as in SDSS J2200-0741) produce secondary minima at 





Table 2 




Best-fit Parameters 


Period 


Frequency 


Amplitude 


(s) 


(mHz) 


{%) 


2008 Jul 27; 






700.7 ± 9.2 


1427 ± 19 


0.55 ± 0.11 


331.7 ± 2.8 


3015 ± 26 


0.40 ± 0.11 


169.30 ± 0.59 


5907 ± 21 


0.50 ± 0.11 


2009 Apr 20: 






337.0 ± 1.4 


2967 ± 12 


0.293 ± 0.080 


169.49 ± 0.41 


5900 ± 14 


0.243 ± 0.081 


2009 Jun 27: 






340.9 ± 1.6 


2933 ± 14 


0.232 ± 0.077 


169.47 ± 0.22 


5900.7 ± 7.8 


0.408 ± 0.077 


2009 Jul 23: 






338.9 ± 1.5 


2950 ± 13 


0.288 ± 0.070 


168.94 ± 0.33 


5919 ± 12 


0.329 ± 0.070 



the location of the fundamental maxima; and smaller 
harmonic amplitudes flatten the fundamental maxima 
(as in SDSS J1426-f5752). We investigate the phase re- 
lationship between the two modes in SDSS J1337-0026 
under the assumption that they are harmonically related 
(if they are not, the shape of the light curve does not 
repeat at the fundamental frequency). To show that this 
assumption is consistent with the data, we compute a 
weighted average of the frequencies on the four nights us- 
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Table 3 

Best-fit Parameters with fi = 2953.6 /^Hz & f2 = 2fi 



Date 


Amplitude^ of fi 


Amplitude^ of f2 


Phase Difference 


(UTC) 


(%) 


(%) 


{^) 


2008 Jul 27 


0.35 ± 0.13 


0.58 ± 0.14 


+ 18 ± 24 


2009 Apr 20 


0.28 ± 0.08 


0.26 ± 0.09 


-12 ± 21 


2009 Jun 27 


0.22 ± 0.08 


0.44 ± 0.08 


-47 ± 23 


2009 Jul 23 


0.29 ± 0.07 


0.33 ± 0.07 


-1 ± 16 



''The amplitudes and their errors have been multiplied by the 
factor given in the text to correct for finite exposure times. 



ing the inverse variances as weights (Ta vloririQQTl ). This 
gives fi = 2957 ± 7 ^Hz and fa = 5906 ± 6/iHz, which 
is consistent with £2 ~ 2fi. 

We refit the data for each night applying this frequency 
constraint. A weighted average of these results gives fi 
= 2953.6 ± 2.7 ^Hz (338.57 ± 0.30 s); the aliasing in 
the combined 2009 amplitude spectrum prevents us from 
confidently determining a more accurate frequency. Con- 
sistent with the data, we assume the frequencies are the 
same on each night and again refit the light curves with 
the frequencies fixed to look for changes in amplitude and 
relative phase. We multiply the best-fit amplitudes and 
their errors by 7rTe2;pf/sin(7rTe3;pf) to correct for the ef- 
fect of a finite exposure time, Texp (jBaldrvlll999D . Table 
[3] lists these results. We report phase difference as the 
number of seconds between the minimum of the harmonic 
and the minimum/maximum of the fundamental and use 
negative values to indicate the harmonic minimum is 
shifted left of the fundamental minimum/maximum. The 
one-sigma errors reported for the phase differences come 
from bootstrap Monte Carlo simulations, and in each 
case the value falls between the sum of the errors for 
the individual phases and the quadrature sum of those 
errors. 

By folding the light curves at the period of the fun- 
damental, we get a picture of these quantitative re- 
sults (Fig. [5]). The pulse shape is like those of SDSS 
J1426-h5752 and SDSS J2200-0741. There is some sug- 
gestion that the phase relationship between the funda- 
mental and harmonic might not be exactly zero, but the 
errors in phase are large making this hard to determine. 
Similarly, there is no statistically significant change in 
amplitude or phase difference among the nights. 

4. DISCUSSION 

The discovery of a fourth variable among the hot 
DQ white dwarf stars means the variable fraction of 
the known hot DQs is at least 36%. This is a high 
percentage compared with the hydrogen- and helium- 
atmosphere white dwarf variables. Thus, as we get to 
know the hot DQ stars better, we find them disinclined 
to be photometrically constant. Further, there is no ob- 
vious observational characteristic that serves as a pre- 
dictor of whether a given hot DQ is variable. Of the 
four reported variables, SDSS J1 426-I-5752 h as a spec- 
tral feature identified with helium (jDufour et al. 2008bJ); 
SDSS J 1337- 0026 has a spectr al feature identified with 
hydrogen (Dufour et al.ll2008ah : two (SDSS J1426-f-5752 
and SDSS J2200-07 41) have broadened spectral lines 
(iDufour et al.ll2008aD pos siblv resulting from a magnetic 
field (jDufour et all2008b[ ): the other two exhibit no such 
spectral distortions. Thus, hydrogen, helium, and mag- 
netic fields do not discourage variability, and none of 
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Figure 2. Light curves for each of the four nights folded at the 
fundamental period. The l-cr errors on the points range from ~ 
0.25% to ~ 0.50%, but the scatter in the data does not repre- 
sent this level of noise since our folding algorithm has the effect of 
smoothing the data. The data are duplicated across two periods. 

them appears necessary to encourage it. Indeed, so far, 
the best predictor of a hot DQ star's variability is that 
it is a hot DQ star. 

Not only do we not know the reason one hot DQ 
varies and another doesn't, we do not kno w why their 
ulse s hapes differ. B a sed on suggestions of IGreen et alj 
20091 ) ■ IDufour et al.l (12009a) predict a connection be- 
tween magnetic field and pulse shape. SDSS J1426-I-5752 
and SDSS J2200-0741 both have relatively large first 
harmonics whose minima coincide with the maxima and 
minima of the fundamental. SDSS J 2 348-09 42 has no ap- 
parent harmonic. IDufour et al.l (2009a) and IGreen et al.l 
(2009) suggest that a magnetic field might account for 
this difference, but SDSS J1337-0026 calls this into ques- 
tion. It has a pulse shape like that of SDSS J1426-f 5752 
and SDSS J2200-0741 but u nlike them its high S/N 
spectrum ( IDufour et al.|[2009b( ) shows no signs of a mag- 
netic field. 

Other connections between the variables' spectro- 
scopic properties and their variable properties are also 
not forthcom i ng. Th e preliminary temperature fits of 
IDufour et al] (l2008al) indicate that SDSS J1337-0026 
is the hottest of the four variables; it is also the 
fastest and smallest in amplitude. However, no obvi- 
ous temperature-period trend emerges when considering 
the other three. Such a relations hip is observed and pre - 
dicted in the ZZ Ceti pulsators (jMukadam et al.l 120061 ). 
Similarly, there is no straight-forward correlation be- 
tween period and amplitude, a relationship also present 
in the ZZ Cetis (pukadam et al..,2006, ). 
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We do know that all four variables have remained 
relatively stable in frequency and amplitude over the 
course of months to a year. The variations observed 
in the present data set — the possible existence of a 
low-frequency peak in the 2008 July amplitude spec- 
trum and the differences in best-fit phase difference and 
amplitude — are not conclusive. 

Resolving these questions will require more observa- 
tions. High S/N photometry of the other hot DQ stars is 
required before any theoretical study can address which 
stars are variable and which are not. SDSS J1337-0026 
was ob served in the original study of IMontgomerv et al.l 
(|2008D but not found to be variable on account of its 
small {'^ 0.3%) amplitude. A mean amplitude spectrum 
noise level of even 0.1% would not permit a convinc- 
ing detection of such small-amplitude variability. We 
also note that the 2009 June data show the harmonic 
at 4.1 times the mean noise in the amplitude spectrum 
while the fundamental is roughly half its size; thus, it 
seems possible that small-amplitude variable hot DQs 
with large harmonics could have fundamentals hidden in 
the noise. 

In addition to high S/N time-series photometry, we 
need high S/N time-series spectroscopy. With higher 
S/N spectroscopy, we might study line profile variations, 
whic h are a diagnostic of non- radial g-mode pulsations 
(e.g.. Ivan Kerkwiik et al.ll2000[ ). th e leading m echanism 
for explaining hot D Q variability (jFontaine e t al. 200l; 
ICorsico et all 120091 ). As with the ZZ Ccti pulsators, 
multi-color photometry will also be an important tool 
for demonstrating that the variability arises prim arily 
from temperature variations (jRobinson et al.|[l982[ ) and 
for determining the spherical harmonics of pulsational 
modes. Unless more than one mode can be detected in 
each star, the prospects for seismology seem dim. How- 
ever, that was also the case for ZZ Ceti stars shortly after 
their discovery, and it is likely that aggressive campaigns 
to detect and observe more hot DQ variables will yield 
fruitful avenues of exploration that are richer and more 
interesting than we could have guessed. 
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